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By generalizing our previous work on the parity symmetry, the partition function of a Yang-Mills 
theory is decomposed into a sum of path integrals each giving the contribution from multiplets of 
states with fixed quantum numbers associated to parity, charge conjugation, translations, rotations 
and central conjugations. Ratios of path integrals and correlation functions can then be computed 
with a multi-level Monte Carlo integration scheme whose numerical cost, at a fixed statistical 
precision and at asymptotically large times, increases power-like with the time extent of the lattice. 
The strategy is implemented for the SU(3) Yang-Mills theory, and a full-fledged computation of 
the mass and multiplicity of the lightest glueball with vacuum quantum numbers is carried out at 
two values of the lattice spacing (0.17 and 0.12 fm). 
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1. Introduction 

Very often the signal-to-noise ratio of correlation functions computed by "standard" Monte 
Carlo techniques decreases exponentially with the time separation of the sources [1,2]. In spec- 
trum computations one has thereby to find compromises between large statistical errors at large 
time-distances and large systematic errors, due to contaminations of excited states, at short time- 
separations. This is not entirely satisfactory from a theoretical point of view since a solid evidence 
that a single state dominates, i.e. a long exponential decay over many orders of magnitude, is 
usually missing. 

The situation is particularly unfavorable in the case of the glueball spectrum, as the variance 
of two-point functions at lai^ge time-distance is dominated by the vacuum [1, 2]. An intuitive way 
to understand the problem starts from the observation that symmetries are usually not preserved 
on a single gauge configuration. No matter what quantum numbers are specified at the source and 
sink, all states are allowed to propagate in the time direction and the expected signal emerges in the 
gauge average as a result of possibly large cancellations. 

In a series of papers [3, 4, 5] (see also [6, 7]) we have proposed and tested a computational 
strategy in which, for each configuration, only states with specified quantum numbers are allowed 
to propagate in the time direction. In such a setup the signal-to-noise problem can be solved 
by introducing a hierarchical integration scheme [8, 9]. We have implemented this "Symmetry 
Constrained Monte Carlo" for computing the mass and the multiplicity of the lightest glueball 
state in the SU(3) Yang-Mills theory. Here we briefly review the basic ingredients entering that 
computation, the results at a rather coarse value of the lattice spacing (a=0.17 fm), and we present 
new numerical data at an additional finer resolution ((3=0.12 fm). 

2. Lattice symmetries 

We aim at computing the ratio of the partition function restricted to a sector, identified by a 
complete set of conserved quantum numbers, over the standard one. To this end, and to fix the 
notation, it is useful to list the symmetry groups of the SU(3) Yang-Mills theory on a finite periodic 
lattice of volume V = T xL^ , where T is its time-extent and L its length in each spatial directions ^ 
We will denote by r^{^i) the matrix associated to the i-th element of the group in the irreducible 
representation jjL. 

• Parity. The group is of order 2, and the two irreducible representations of dimension 1 are 

r(±)(=^i) = ±rW(^2) = i- 

• Charge Conjugation. Again this group is of order 2, and the two irreducible representations 
of dimension 1 are r^^\^i) = ±r(^'(^2) = 1- 

• Translations. The group of translations is a direct product of three Abelian groups, one 
for each space direction. Its elements are labeled by a three dimensional vector of integers 
m = {mi, 1112, m^), with m, = 0, . . . ,L — 1, where each component identifies the elements of 
the group in the corresponding direction. Since it is Abelian, each element forms its own 
class and there are non-equivalent irreducible representations of dimension 1 

r(p)(=^„)=f''>™, (2.1) 

^Dimensionful quantities are always expressed in units of a unless explicitly specified. 
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which are labeled by momentum vectors p= — [?ii,«2,«3],with?i,=0,...,L— 1. 

• Rotations. The octahedral group is of order 24. Its elements are listed in Appendix B of [5], 
where explicit expressions for the Y matrices and the table of characters can be found for the 
5 non-equivalent irreducible representations. Those are two singlets Ai and A2, one doublet 
E and two triplets Ti and T2. 

• Central Charge Conjugations. This symmetry is strictly related to the choice of periodic 
boundary conditions, and it disappears in the infinite volume limit [10]. The group is a direct 
product of three Z3, one for each spatial direction. It is of order 27, and its elements are 
labeled by a three dimensional vector of integers v = ( Vi , V2 , V3 ) , with v,- = , 1 , 2 , where each 
component labels the elements of the Abelian group in the corresponding direction. Since 
each element forms its own class, there are 27 non-equivalent irreducible representations of 
dimension one 

^W(^v)=^'''■^ (2.2) 
27r 

which are labeled by the electric flux vectors [^1,^2, ^3], with ej = 0, 1,2. 

3. Symmetry constrained Monte Carlo 

For a discrete group of order g, the projector onto states which transform as an irreducible 
representation /i can be defined as (for unexplained notation see Ref. [5]) 

where is the dimension of the iiTcducible representation, X^^\^i) is the character of the ith 
group element in that representation, and f (^;) is the representation of the element in the Hilbert 
space. The corresponding symmetry-constrained partition function can then be expressed as 

= Xi.|t^p/i} ^ (3 2) 

where T is the transfer matrix among gauge-invariant states. By inserting (3.1) in the equation 
above and by choosing the "coordinate" basis to express the trace, it is clear that Z^^^\T) can be 
written as a linear combination of partition functions of g different systems with twisted boundary 
conditions in the time direction. The latter are chosen so that the state at time T is related to the 
one at time by a group transformation. As shown in detail in Ref. [4], the ratio Z^'^' (T) /Z{T) can 
be factorized as a product of similar ratios associated to thick time-slices of temporal extension d 
and fixed boundary conditions. Explicitly, we numerically compute 



z(A')fr) 1 



z(r) z(r) 



D«76'-^[^]p';f][r,0], (3.3) 



where Pf,f^[r,0] is a product ofm = T /d factors, and depends on the values of the spatial links on 
the boundaries of the thick time-slices only. Its definition is given in Refs. [4, 5] . Finally, given the 
locality of the gauge theory, such a factorized quantity can be very accurately estimated through 
a generalization of the hierarchical integration scheme proposed in Refs. [8, 9], which in this case 
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Table 1: Simulation parameters: A^conf is the number of configurations, A^iev is the number of levels and d is 
the thickness of the thick time-slice used for the various levels. 

removes completely the exponential signal-to-noise problem [3, 4, 5]. Although we focus here on 
partition functions, the same approach can be applied to the computation of correlation functions. 
Matrix elements of operators among glueball states can again be obtained avoiding the exponential 
signal degradation [5]. 

To determine the mass of the lightest glueball state with vacuum quantum numbers, we are 
interested in computing the ratio z(*"*''P'''^=+)(r) /Z{T), which we will shorten asZ(P'+^(r) /Z{T). 
The projection onto non-zero momentum is needed to get rid of the contribution from the vacuum, 
which would otherwise dominate the variance and cause an exponential degradation of the signal 
from the glueball. By applying blindly the analysis in this section, we would need to calculate 
the thick time-slice ratios for each of the x 27 x 2 boundary conditions. This would make the 
approach extremely expensive from the computational point of view, and it would give unnecessary 
information if one is interested in the low momentum states only. It is possible, however, to have 
still an exact numerical algorithm by implementing the projectors on and ^ exactly, while 
treating those on {py,P:) = (0,0) and e = stochastically by extracting n random transformations 
out of the associated x 27 ones and averaging over them, see Ref. [5]. As we fix singlet quantum 
numbers for py , p^ and e, a stochastic treatment seems a priori justified, fluctuations are proportional 
to the exponentially suppressed higher momentum components, or to torelon contributions, which 
are expected to have higher energies on these volumes. The results in the next section have been 
obtained by using values of n between 9 and 64, a choice justified a posteriori by the moderate 
statistical errors obtained. 

4. Results 

We have simulated the SU(3) gauge theory discretized on the lattice by the Wilson action at 
j3 = 6/gQ = 5.7, and j8 = 5.85 which correspond to a spacing of 0.17 fm and 0.12 fm respec- 
tively [11]. The spatial lengths are 1.4 and 1.7 fm, while time extends up to 2 fm. The simulation 
parameters and the results are summarized in Tables 1 and 2. 

The primary quantity that we have computed is the ratio Z(P'+) /Z with 772.3 = 0, pi = {2n/L)n[ 
and ni = 1,2. The results for the A series ai^e shown in the left plot of Figure 1. We fit them by 
assuming that a single state contributes, i.e. using an ansatz of the form ln(Z(P'+)(r)/Z(r)) = 
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1.0(3) -10-3 


1.15(5) 
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0.94(25) • 10-4 


0.92(25) • 10-4 


1.55(5) 


Ci 
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1.5(3) -10-4 


1.5(4) -10-4 


0.883(24) 
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2.4(22) -10-5 


2.4(22) • 10-^ 





Table 2: Results for ratios of partition functions with momenta p ~ [2;r«i/L.0,0]. The effective energy 
^eff is defined as in Eq. (4.1). 




Figure 1: Left panel: ratio of partition functions z'P +^/z''''+' with momenta p = [27r/L,0,0] for the A 
lattices. Right panel: the corresponding effective energy as defined in Eq. (4.1). The band is our best 
estimate, i.e. the one extracted from the lattice with the longest time-extension. 



A — BT , such that A yields the logarithm of the multiplicity of the state and B its effective, finite 
momentum, energy E^^^\ At lai^ge time separations the lightest glueball with vacuum quantum 
numbers and momentum p is expected to dominate. The four points at the largest values of T 
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(r/fl; = 5,6,8,12)in Figure 1 ar^e well described by our ansatz, and the fit results for the multiplicity 
are well consistent with a value of 1 excluding 2 and 3 by three and four standard deviations 
respectively. We therefore define 

^eff — J ^" 

for which the results ai^e summarized in Table 2 and shown on the right-hand plot of Figure 1. 

The lattice B3 serves the purpose of assessing finite size effects. It has the same lattice spacing 
of the A series but a linear- extension of L = 10. The results for «i = 1,2 are reported in Table 2, 
and ai^e plotted as a function of the momentum squared in the left plot of Figure 2. The dashed 
line is a linear interpolation of the two black points (circles) of the lattice B3, while the red point 
(square) is our best result for the A series. It is rather clear that, within our statistical precision, 
finite volume effects are not visible in our data. It is also interesting to notice that, even if the values 
of the momenta are rather large, the continuum dispersion relation is well reproduced within our 
statistical errors. We extract the mass of the lightest glueball from the expression 

which, using the T /a = \2 result and in units of the lattice spacing, gives 

M+= 0.935 ±0.042, [jS = 5.7] , (4.3) 

in good agreement with the estimate in Ref. [12] obtained at the same lattice spacing with the same 
discretization but within the standard approach. Finally, the lattice Ci is matched to B3 in volume 
but with a larger time extension (corresponding to T/a = 7 at j3 =5.7) and more importantly with 
a finer lattice resolution, namely a = 0.12 fm. By making use again of Eq. (4.1) to extract E^^^^\ 
we get from the continuum dispersion relation 

M+ =0.760 ±0.028 , [/3=5.85], (4.4) 

in units of the lattice spacing. The results, measured in units of the scale ro [11], ai^e collected in 
Fig. 2, where they are plotted as a function of {a/rof- as the leading discretization effects should 
be quadratic in the lattice spacing. 

5. Conclusions 

We have discussed how the relative contributions to the partition function, due to states carry- 
ing a given set of quantum numbers associated with the exact symmetries of a field theory, can be 
expressed by ratios of path integrals with different boundary conditions in the time direction. From 
an algorithmic point of view, the composition properties of the projectors can be exploited to imple- 
ment a hierarchical multi-level integration procedure which solves the problem of the exponential 
(in time) degradation of the signal-to-noise ratio. 

Within this approach we have performed a precise lattice computation of the mass of the 
lightest glueball with vacuum quantum numbers in the SU(3) Yang-Mills theory at two values of 
the lattice spacing corresponding to 0.17 and 0.12 fm. The algorithm works as expected, and we 
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E" = 0.9G(15) + 0.51(15)* 



0.2 0.4 06 OR 



3.6 - 
3.4 - 
3.2 - 
3 - 
2.8- 
2.6 - 
2.4 - 
2.2 - 



0.02 0.04 0.06 O.Oi 



0.1 0.12 0.14 



Figure 2: Left panel: the effective energy squared from the A5 (red square) and the B3 (black circles) lattices. 
Right panel: results for roM+ at a ~ 0.17 fm and a ~ 0.12 fm. 



have been able to follow the exponential decay up to separations of 2 fm, while keeping the error 
on the effective mass approximatively constant as a function of time. 

Cutoff effects appear to be rather large and tend to significantly decrease the estimate of the 
glueball mass at finite lattice spacing, as expected and also observed in previous lattice computa- 
tions [13, 14]. We are now in the process of repeating the calculation presented at a finer lattice 
resolution of 0.1 fm, which should eventually allow us to properly assess the magnitude of dis- 
cretization effects. 

The simulations were performed at CILEA, at the Swiss National Supercomputing Centre 
(CSCS) and at the Jiilich Supercomputing Centre (JSC). We thankfully acknowledge the computer 
resources and technical support provided by all these institutions and their technical staff. 
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